The triple alpha reaction rate and the 2 + resonances in 12 C 



R. de Diego and E. Garrido 
Institute) de Estructura de la Materia, CSIC, Serrano 123, E-28006 Madrid, Spain 



o 

> 
O 

(N 



O 
=5 

S3 



> 
O 

in 
in 
o 



X 



D.V. Fedorov and A.S. Jensen 
Department of Physics and Astronomy, Aarhus University, DK-8000 Aarhus C, Denmark 

(Dated: January 12, 2011) 

The triple alpha rate is obtained from the three-body bound and continuum states computed in 
a large box. The results from this genuine full three-body calculation are compared with standard 
reference rates obtained by two sequential two-body processes. The fairly good agreement relies on 
two different assumptions about the lowest 2 + resonance energy. With the same 2 + energy the rates 
from the full three-body calculation are smaller than those of the standard reference. We discuss the 
rate dependence on the experimentally unknown 2 + energy. Substantial deviations from previous 
results appear for temperatures above 3 GK. 

PACS numbers: 25.10.+S, 25.40.Lw, 26.30.Hj 



Introduction. The triple alpha process is the key re- 
action that permits to bridge the A — 5 and ^4 = 8 gaps, 
opening the door to the production of 12 C in the core of 
the stars in the red giant phase 

The properties of the a-a-a continuum states are cru- 
cial for the reaction rates, which determine the abun- 
dance of 12 C in the Universe. As a clear example, this 
is what led Fred Hoyle to predict the existence of a + 
resonance just above the triple a threshold Q • It is not 
possible to explain the observed abundance of 12 C with- 
out this so called Hoyle state, which was confirmed ex- 
perimentally 0] to be at an energy of 0.38 MeV above the 
three-body threshold. The Hoyle state enhances by itself 
the reaction rate by about two orders of magnitude at low 
temperatures (below 2 GK), where the rate is dominated 
by the electric quadrupole transition from continuum 3a 
+ states to the excited bound 2 + state in 12 C. 

In type-II supernova explosions dense and hot envi- 
ronments are created. This is the so called hot bubble, 
consisting of alpha-particles and neutrons, with rather 
uncertain but relatively high temperature reaching sev- 
eral GK [H . These environments are the suggested place 
for the rapid-neutron process Although the impor- 
tance of the triple alpha contribution as compared to 
other reactions in this temperature range (T > 2 GK) 
has been questioned @, it has also been suggested that 
changes in this rate imply changes in estimates of the rel- 
ative amounts of elements formed during the supernovae 
explosion, and therefore in estimates of the rate at which 
heavy elements are distributed through the Universe Q- 
In [8| the estimated reaction rate for the triple alpha re- 
action in this temperature range was found to be about 
an order of magnitude lower than the one given in 
which would lead to a mass fraction of 56 Ni about two to 
three times smaller. 

At temperatures above 3 GK, the reaction rate for 
the triple alpha reaction is dominated by the electric 
quadrupole 2 + — > Of transition [9] (we follow the no- 
tation where J* represents the i th state in the spectrun 
with angular momentum J and parity n. In the partic- 



ular case of 12 C the states Of and 2^ are the only ones 
bound). The reaction rate in this region is therefore to a 
large extent determined by the properties of the 2 + con- 
tinuum states, and in particular by the possible existence 
of a 2^~ resonance at a relatively low energy. The proper- 
ties of such a state is still an open problem. The reaction 
rates given in [§, [l(| were computed assumingthat 12 C 
has a 2^~ resonance at 1.75 MeV. However, in Q no evi- 
dence was found concerning the existence of such 2 + res- 
onance. On the contrary, in (Tl| an experimental energy 
of about 2.3 MeV is derived. This value is substantially 
higher than the one used in [j| [l(| , and it could substan- 
tially modify their computed results. Also, in [lj|], where 
a three-body calculation is employed to obtain the reso- 
nances as poles of the S'-matrix, a 2^~ resonance has been 
found at an energy of 1.38 MeV. However, this energy is 
found with the same (adjustable) three-body interaction 
that reproduces the energy of the bound 2^ state. 

The role of the 2^~ resonance in 12 C in the triple alpha 
reaction rate requires detailed investigation. We shall 
employ the same three-body method as in 13], where no 
assumption is made about the capture mechanism (se- 
quential or direct), and in contrast to the method de- 
scribed in Q, where a sequential capture process is as- 
sumed. The standard procedure and the almost canoni- 
cal results of 9] is furthermore tested by comparison. 

Full three-body formulation. The first of the methods 
used in this work is the three-body calculation described 
in (T^ |. Let us consider the radiative capture process 
a + b + c — > A + 7, where A is a bound system made 
of particles a, b, and c with separation energy Ba ■ The 
corresponding reaction rate R a bc(E) is given by [H, [H| 
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where E = E 1 + Ba is the initial three-body kinetic 
energy, E y is the photon energy, cx 7 (i? 7 ) is the photo dis- 
sociation cross section of the A nucleus, c is the velocity 
of light, gi is the spin degeneracy of states of particle 
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i = a,b, c, A, and fx x and \i y are the reduced masses of 
the systems related to the Jacobi coordinates, (x, y), for 
the three-body system [la ]. 

The photo dissociation cross section for the inverse pro- 
cess A + 7 — > a + b + c can be expanded into electric and 
magnetic multipoles. In particular, the electric multipole 
contribution of order A has the form Fiji : 
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where the strength function B is 
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B(EX,n Q J -> nJ) = ^ |(nJM|O^|n J M }| 2 , (3) 

where Jo, J and Mo, M are the total angular momenta 
and their projections of the initial and final states, and 
all other quantum numbers are collected into no and n. 
The electric multipole operator is given by: 
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where i runs over the three particles of charges Zj , and 
where we neglect contributions from intrinsic transitions 
within each of the three constituents [l6| . 

Finally the energy averaged reaction rate is obtained 
as a function of the temperature by using the Maxwell- 
Boltzmann distribution as weighting function. For three 
alpha particles we obtain [l3j : 
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where k B is the Boltzmann constant. 

The strength function B is computed by genuine 
three-body calculations of both the bound final state, 
noJoMo), and the continuum initial states, \nJM). We 
use the hyperspherical adiabatic expansion method de- 
scribed in [15j . The a-a interaction is given in [l7j . The 
basic procedure is computation of three-body states of 
given angular momentum and parity confined by box 
boundary conditions [l8j]. In this way the continuum 
spectrum is discretized. The strength functions are then 
obtained for each discrete continuum state according to 
Eq.(|3|). The distribution dB/dE is built by use of the 
finite energy interval approximation, where the energy 
range is divided into bins, and all the discrete values of B 
falling into a given bin are added. Afterwards the points 
are connected by spline operations and the expressions 
© and ([5]) arc computed. 

Sequential process. The NACRE results given in Q 
are usually taken as the reference for the reaction rates 
in the triple alpha process. In this work the process is 
assumed to proceed in a sequentially two-step process. 



In the first step one a-particle captures another one to 
produce 8 Be in the ground + resonant state. In the sec- 
ond step, 8 Be is able to capture (before decaying) a third 
a-particle, populate a 12 C resonance, and then decay by 
photo emission to one of the bound states of 12 C. The re- 
action rate for such a two-step process is given by the rate 
for the capture of an a-particle by 8 Be {{RcfiaJJE" , E')}) 
weighted with the rate for formation of 8 Be [9|: 
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where E" and E' are the relative energy between the 
two a-particles forming 8 Be and the energy of the third 
a-particle relative to the center of mass of the first two, 
respectively. The temperature is T and fi aa is the re- 
duced mass of the two-alpha system. The elastic a — a 
cross section is given by: 

2ir T a ( 8 Be,E") 2 
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where k 2 = 2/j, aa E" /h 2 , E r is the 8 Be resonance energy, 
and the width r Q ( 8 Be, E") has the form [§]: 
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where T a is the width of the resonance, Pi is the pene- 
tration factor, and t is the relative orbital angular mo- 
mentum between the two a-particles. 

In our case, the experimental energy of the + reso- 
nance in 8 Be is E r = 0.09189 MeV above threshold, with 
a width of T a = 6.8 ± 1.7 ■ 10~ 6 MeV. For such a narrow 
resonance we have that r Q ( 8 Be,_E") « T Q , and for our 
purpose we can safely replace Eq.Q by: 
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Finally, from Ref.[9j we have: 
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where 
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TABLE I: Values of the parameters entering in Eqs. (|12[l to 
(|15|) taken from Ref.[9j. They are all given in MeV. 



where is the 12 C resonance energy with angular mo- 
mentum J, and 
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The parameters needed to compute Eq. (fl2j) . and there- 
fore (fTT|) and (fT0|) . are taken from [9( and they are col- 
lected in table HI 

Given a transition between 12 C continuum states with 
angular momentum J and a bound 12 C state with bind- 
ing energy Bj, it is now possible to extract the expression 
for dB/dE effectively used in Ref.@]. This can be made 
by inserting Eq.Ql]) into (JTOJl, and Eq.® with A=2 into 
Eq.©. By comparison of those expressions we extract: 
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where E = E r + E' is the three-body kinetic energy and 
a is the fine structure constant. 

Comparing the methods. In [l3| the full three-body 
method was used to compute the reaction rate for the 
triple a process. The a-ct interaction in (l7| was used. 
This leads to the 12 C spectrum given in [12j, which results 
in good agreement with the experimental data for about 
13 different states. In particular, concerning the 2 + reso- 
nances, the computed 12 C spectrum includes a low lying 
one at 1.38 MeV (above the three-body threshold) with 
a computed width of 0.13 MeV, and a second one at 4.4 
MeV with a width of about 1 MeV. Experimentally only a 
resonance at 3.88 MeV is fully established. In Ref.[12| the 
resonances are obtained as poles of the iS-matrix by use of 
the complex scaling method [l9| . With this method the 
resonances behave asymptotically as bound states with 
complex energy, which automatically determines the res- 
onance energy and the corresponding width. The res- 
onance width is then not a parameter but an output of 
the calculation. The three-body calculation, as described 
in [l]|, gives rise to the triple a reaction rate shown by 
the thick solid line in FigQ] The contributions coming 
from the + — > 2+ and 2 + — > Of transitions are shown 
by the corresponding thin solid lines. The agreement be- 
tween the total reaction rate and the result given in Q 
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FIG. 1: Reaction rate for the triple a process with the full 
three-body calculation (thick solid line) and the sequential 
approach (thick dashed line) as described in the text. The 
corresponding contributions from the + — > 2+ and 2 + — > Of 
transitions are given by the thin curves (2~f and 0^ refer to 
the lowest 2 + and + states in the 12 C spectrum, which are 
the only ones bound). The thin dotted line is the contribu- 
tion from the 2 + — > 2f in the full three-body calculation. 
The open circles correspond the rate given in The thick 
dotted curve is a sequential calculation but using the same 
energy and width for the 2j resonance as in the three-body 
calculation. 



(open circles in the figure) is reasonably good. For com- 
pleteness, we show in the figure the contribution from the 
2 + — > 2~i transition (thin dotted curve). This contribu- 
tion is very small and could actually be neglected. 

When the sequential approach is used, and the reaction 
rate is obtained from Eq. (fT0|) . we get the thick dashed 
curve in the figure. The sequential method described in 
the present work and in are identical except for the 
approximation in Eq.©. The parameters used in the cal- 
culations (table []| are also the same in both cases. The 
good agreement between the thick dashed curve and the 
circles in FigQ] shows that the approximation in Eq.© 
is accurate. Only a small difference is found for high 
temperatures since the contribution from the 3j~ reso- 
nance in 12 C is included in [9j, but completely neglected 
in our calculation. Again, the corresponding thin dashed 
lines show the contributions to the total rate from the 
0+ — > 2^ and 2+ — > Of transitions. 

As seen in FigQ] the contribution from the + — > 2f 
transition is very similar in the full three-body calcula- 
tion and the sequential approximation. This is due to the 
narrow 0^" Hoyle three-body state in 12 C which heavily 
dominates the full calculation. It has a large strength 
corresponding to strong population and subsequent de- 
cay into the bound 2f state plus a photon. This process 
is precisely as assumed in the sequential description, and 
in fact the reason for its success. 

However, the lowest 2^~ resonance in 12 C (both the one 
obtained in the three-body calculation and the one used 
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in Q) is rather wide, which means that the continuum 
non-resonant three-body 2 + states can contribute signif- 
icantly to the 2 + — > Of transition. In fact, as seen in the 
figure, the three-body and the sequential calculation dif- 
fer for this contribution, especially at low temperatures. 
Due to this discrepancy the two computed total reaction 
rates do not fully agree for temperatures from about 2 
to 5 GK. When the temperature increases the thin solid 
and thin dashed curves in FigfT]become closer and closer, 
in such a way that for temperatures higher then 5 GK 
the total reaction rate is similar in both calculations. 
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FIG. 2: Upper part: the dB/dE function for the 2+ -> 0+ 
transition for the sequential approach as given in Eq. (|16[) 
(dashed curve) and the three-body calculation (solid curve). 
Lower part: Integrands of Eq.© (thick curves) and Eq. (|10[) 
(thin curves) in arbitrary units as a function of the kinetic 
energy of the three q's for T=3, 5, 7, and 9 GK. The thick 
lines and thin lines correspond to the full three-body and the 
sequential calculations, respectively. 



The fact that both calculations produce basically the 
same result is encouraging. However, this good agree- 
ment is apparently inconsistent, since it is achieved in 
spite of different energies of the lowest 2\ resonance in 
12 C. While in the full three-body calculation this energy 
is 1.38 MeV, the sequential result is computed for an 
energy of 1.75 MeV. In fact, when the sequential calcula- 
tion is performed using the same energy (1.38 MeV) and 
width (0.13 MeV) as in the full three-body calculation we 
obtain the rate given by the thick dotted curve in FigJT] 
As we can see, for the same parameters in the low-lying 
2^~ resonance in 12 C, the sequential calculation gives a 
rate that, for temperatures beyond 2 GK, is clearly above 
the one obtained in the three-body calculation. 

For a better understanding it is useful to compare the 
dB/dE function for the 2 + — > 0^ transition obtained 
in the full three-body calculation and in the sequen- 



tial approximation from Eq.([T6]). This comparison is 
shown in the upper part of FigJ5]by the solid and dashed 
curves, respectively. The dB/dE function in the sequen- 
tial case is essentially a lorentzian centered around the 
1.75 MeV corresponding to the 2 J resonance energy. In 
the full three-body calculation the peak, as expected, is 
located at about 1.4 MeV. Furthermore, a shoulder ap- 
pears for higher energies arising from the contribution 
from continuum non-resonant states and the additional 
2 + resonances obtained in the three-body calculation. 
These contributions are not included in the description 
in Ref.|H. 

However, due to the exponentials in Eq.© or (fTTj) . the 
high energy region of the transition strength contributes 
very little to the reaction rate. This is shown in the lower 
part of the figure, where we plot the integrands entering 
in the calculation of the reaction rate (Eqs.([5]) and (TTO)) ) 
as a function of the three-body kinetic energy for T=3, 
5, 7, and 9 GK. The thick and thin curves correspond 
to the full three-body and the sequential calculations, 
respectively. As expected, these integrands die rather 
fast with the energy, and the larger the temperature the 
larger the energy range contributing to the reaction rate. 

As seen in the upper part of the figure, the full three- 
body calculation accumulates part of the strength at 
smaller energies than the sequential calculation. This 
is first of all due to the different resonance energies of 
1.38 MeV in the three-body and 1.75 MeV in the sequen- 
tial calculation. For this reason, for low temperatures, i.e. 
for T < 3 GK, where the integrand is vanishingly small 
for E > 2 MeV (see the solid curves in the lower part of 
the figure), the three-body calculation provides a higher 
reaction rate than in the sequential approximation. The 
integral of the thick solid line in the lower part of Fig|2] 
is almost twice the one of the thin solid curve. When the 
temperature increases, the energy at which the integrand 
begins to vanish also increases, and the sequential reac- 
tion rate adds more strength compared to the three-body 
case. As a consequence both reaction rates approach each 
other, and they even cross at some point. The integral of 
the thin dot-dashed curve in the lower part of Fig(5] (se- 
quential case for T=9 GK) is about 25% bigger than the 
one of the thick dot-dashed curve (three-body case for 
T=9 GK). Even for T = 9 GK the integrand only allows 
contribution up to energies around 5 MeV, which means 
that the behavior of the dB/dE function for E > 5 MeV 
is rather unimportant. 

Thus, from FigJ5]we can conclude that the full three- 
body calculation shifts part of the strength to large ener- 
gies, which contribute very little to the reaction rate. For 
this reason, for the same low lying 2^" resonance energy 
in 12 C, the full three-body calculation provides a smaller 
rate than the sequential calculation. Or, in other words, 
for the three-body and sequential calculations to provide 
the same reaction rate, the 2^~ resonance energy has to 
be smaller in the three-body case, in order to compensate 
the non-contributing part of the strength at large ener- 
gies (for instance 1.38 MeV in the three-body calculation 
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and 1.75 MeV in the sequential approach). 

The 2t resonance energy and the reaction rate. As 
seen in Fig[TJ for temperatures below about 3 GK the 
reaction rate of the triple a reaction is dominated by the 
+ — > 2f transition, and the precise energy of the lowest 
2^~ resonance in 12 C does not play any role. Only for tem- 
peratures beyond 4 GK the energy of the 2^~ resonance 
has a sizable effect. 
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FIG. 4: Reaction rate for the three-body calculation when 
the 2j resonance is placed at 1.38 MeV (thick solid line), 
when the resonance is removed from the calculation (thick 
dashed line), when the full contribution from the 2 + — \ 0^ 
transition is excluded (thin solid line), and when the Hoyle 
resonance is removed (dotted line). The thin dashed curve is 
the calculation in the sequential case when the energy of the 
2 J resonance is 4.0 MeV. The open circles are taken from [{J. 



FIG. 3: Reaction rate in the sequential case for different ener- 
gies of the lowest 2% resonance in 12 C. The energy increases 
from the upper curve to the lower one from 1 MeV up to 
5 MeV. The dotted curve is the calculation where the con- 
tribution from the 2 + — >• 0^ transition has been completely 
removed. 

FigfS] shows the total sequential reaction rate for dif- 
ferent energies of the 2^~ resonance in 12 C. The rate de- 
creases progressively when the resonance is placed at 
a higher and higher energy. This happens because the 
higher the energy of the 2^~ resonance, the more strength 
located at high three-body energies, where the integrand 
is small due to the exponential in Eq.© or (|TT]) (see 
the lower part of Fig 12]). In fact, for a resonance energy 
of 5 MeV the computed rate is quite similar to the one 
obtained when the resonance is removed. Only for tem- 
peratures beyond around 7 GK some difference appears. 

For the same reason, the same behavior is found for the 
three-body calculation. An increase of the 2^~ resonance 
energy decreases the reaction rate. However, full removal 
of the 2^~ resonance in 12 C gives a rate that differs quite 
a lot from the one obtained when the full contribution 
from the 2 + — > 0^ transition is removed. This can be 
seen in Fig 01 where the thick solid line is the same cal- 
culation as in Fig[T](with the 2 J resonance at 1.38 MeV), 
the thick dashed line is the calculation after removal of 
this 2\ resonance, and the thin solid line is the result 
when the full 2 + — > 0^ transition is excluded. The dif- 
ference between the thin solid line and the thick dashed 
line in the figure shows the importance of including the 
additional continuum 2 + states in 12 C. 

It is important to keep in mind that the three-body 
calculation, together with the 2^~ resonance at 1.38 MeV, 



includes as well all the other states obtained in [12|, in 
particular the 2 3" resonance around 4 MeV. For this rea- 
son, the thick dashed line in FigUJ where the resonance 
at 1.38 MeV has been suppressed, still contains the con- 
tribution of the one at about 4 MeV, which is the one 
experimentally known. Then the thick dashed line in the 
figure gives a lower limit to the triple a reaction rate, 
since this is what we would get in case that the 2 + reso- 
nance in 12 C around 1.5-2.5 MeV is shown not to exist. 

In this connection it is interesting to compare the thick 
dashed line in the figure with the thin dashed line, which 
represents the reaction rate obtained in the sequential ap- 
proach when the lowest 2 J resonance in 12 C is placed at 
4.0 MeV. This is very close to the energy of the lowest 2 j" 
resonance fully confirmed experimentally, and therefore 
this would be what provided by the sequential calcula- 
tion in case that the 2 + resonance in 12 C around 1.5-2.5 
MeV does not exist. The difference between both dashed 
curves shows the effect of the non-resonant 2 + states in 
this case. 

For illustration we also show in Fig0]the reaction rate 
that we obtain when the Hoyle resonance is removed from 
the calculation (dotted curve in the figure). As we can 
see, this state is responsible for an increase of the reaction 
of about two orders of magnitude for temperatures below 
3 GK. 

Thus, from FigU we can conclude that in case that 
the low-lying l\ resonance is confirmed not to exist 
(such that the lowest one is the already known resonance 
close to 4 MeV above threshold) the sequential approach 
would underestimate the triple- a reaction rate. The non- 
resonant 2 + states in 12 C are enough to increase signifi- 
cantly the rate at high temperatures. 
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Summary and conclusions. We investigate the reac- 
tion rate for the triple a reaction. We used the full three- 
body calculation described in [l3[ and the sequential de- 
scription of [9] . For the same set of three-body resonances 
we find that the full three-body calculation gives smaller 
rates than the sequential approximation for temperatures 
above ~ 3 GK. 

We focus on the importance of the experimentally un- 
certain lowest 2~£ resonance in 12 C which to a large ex- 
tent determines the reaction rate for temperatures above 
3 GK. Different theoretical calculations predict a res- 
onance in 12 C at around 1.0-2.5 MeV. We find that when 
the full and sequential calculations are performed tak- 
ing a 2^~ resonance energy of 1.38 MeV and 1.75 MeV, 
respectively, the computed total reaction rates are quite 
similar for the whole range temperature. The only excep- 
tion is for temperatures ranging between about 2 and 5 
GK, where the full calculation gives a rate slightly above 
the one obtained with the sequential calculation. The 
reason for this behavior is that part of the strength in 
the three-body case is moved to higher energies due to 
the higher 2 + resonances and non-resonant continuum 
2 + states. The high energy part of this strength con- 
tributes very little, and the missing low-energy strength 
has to be compensated by a lower-lying 2 + resonance. 
For the same energy and width of the 2 resonance, the 
full three-body calculation contributes less than the se- 
quential approximation. 



If the lowest 2% resonance in 12 C by chance should be 
at about 4 MeV, where the lowest established 2 + state 
is located, then the reaction rate is higher for the full 
three-body calculation than obtained from the sequen- 
tial approximation. This is due to contributions from the 
non-resonant continuum states in the three-body calcula- 
tion and the insignificant contribution from a high-lying 
resonance in the sequential calculation. The numerical 
result can be taken as a lower limit to the reaction rate, 
which for high temperatures can be up to about one order 
of magnitude above the result obtained in the sequential 
picture. 

Summarizing, a detailed description of the triple a re- 
action rate for temperatures beyond 3 GK requires a 
careful treatment of the 2 + — > 0^ transition. The non- 
resonant continuum 2 + states, not included in the se- 
quential description, are very important. For a given 
energy of the lowest 2j~ resonance in 12 C below about 
3 MeV, a sequential description of the process overesti- 
mates the reaction rate. For a larger resonance energy 
the rate of three-body calculation exceeds that of the 
sequential approximation due to the contribution of the 
non-resonant continuum states. 
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